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The ground-state phase diagram of Heisenberg spin-1/2 system on a two-leg ladder with rung 
alternation is studied by combining analytical approaches with numerical simulations. For the case 
of ferromagnetic leg exchanges a unique ferrimagnetic ground state emerges, whereas for the case of 
antiferromagnetic leg exchanges several different ground states are stabilized depending on the ratio 
between exchanges along legs and rungs. For the more general case of a honeycomb-ladder model 
for the case of ferromagnetic leg exchanges besides usual rung-singlet and saturated ferromagnetic 
states we obtain a ferrimagnetic Luttinger liquid phase with both linear and quadratic low energy 
dispersions and ground state magnetization continuously changing with system parameters. For 
the case of antiferromagnetic exchanges along legs, different dimerized states including states with 
additional topological order are suggested to be realized. 


INTRODUCTION 



S2,2j-n 


Spin-l/2 Heisenberg two-leg ladder systems have at¬ 
tracted a great deal of interest both from experiment 
and from theory [ij. Ladders with antiferromagnetic ex¬ 
change along run^ and antiferromagnetic as well 
as ferromagnetic ]1, exchanges along legs have been 
realized experimentally. 

For the case of a ladder with ferromagnetic legs and a 
uniform inter-leg (rung) exchange, so called rung-singlet 
or saturated ferromagnetic phases are realized depend¬ 
ing whether the inter-leg coupling is antiferromagnetic 
or ferromagnetic. 

For the case of a ladder with antiferromagnetic legs and 
a uniform rung coupling it has been established that for 
antiferromagnetic inter-leg coupling a rung-singlet phase 
is realized, whereas for ferromagnetic rung coupling a 
Haldane phase is stabilized. Both phases are stable both 
in weak rung-coupling and in strong rung-coupling lim¬ 
its. Weak rung-coupling limit is a proper limit for effec¬ 
tive field theory bosonization analyses Q, where in the 
case of uniform rung exchanges, in the lowest (first) or¬ 
der of the inter-chain coupling the relevant operators (in 
the renormalization group sense) are present that at low 
energies drive system towards the strong-coupling fixed 
points of the rung-singlet and Haldane states for positive 
and negative rung exchanges respectively. 

Stability of a unique ground state from arbitrary weak 
(non-zero) up to arbitrary strong inter-leg exchanges, is 
due to the fact that the ladder system is non-frustrated. 
The question that we are going to address in our work 
is what happens when the ladder system is frustrated 
by rung exchange alternating in sign from rung to rung. 
Frustration in this case, for both signs of exchanges along 
the ladder legs, will be caused by the presence of an odd 
number of antiferromagnetic exchanges in the elementary 
closed path that is a ladder plaquette in our case. This 
problem has not been addressed before and we will try 
to fill this gap in the following. 
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FIG. 1. Geometry of the two-leg ladder with alternated 
rung exchanges. Antiferromagnetic couplings along are cho¬ 
sen along the even rungs, J± > 0. 

We will use different complementary analytical ap¬ 
proaches: strong-rung coupling expansion for strongly 
coupled legs and bosonization for weakly coupled antifer¬ 
romagnetic legs. To cover the intermediate regimes we 
will use numerical techniques. We will as well consider 
the more generalized case of a rung-alternated model 
where we will relax the constraint of equal absolute value 
of exchanges along the even and odd rungs. 

MODEL OF FRUSTRATED SPIN LADDER 

In this work we study a Heisenberg spin-1/2 model 
defined on a two leg ladder with L rungs and with alter¬ 
nating rung exchanges, depicted in Fig. [U 

L L 

h = jJ2 + Ji- E(-i)"Sb.S2,,, (1) 

3 = 1 , 1 - 3 = 1 

where Sij are spin-^ operators acting on spins on the 
j-th rung of the I = 1,2 leg. For definiteness we will put 
Jj_ > 0, the case of Jj_ < 0 will be recovered by one-site 
translation of the ladder along the legs. 

The anisotropic XY case of Eq. m was studied 
recently in the context of single-component hard-core 
bosons on a two-1^ ladder at half-filling with a flux tt 
per plaquette [1, [J. It was shown that depending on 
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the ratio of strengths of exchanges along ladder legs and 
rungs there are two different ground-state phases. For 
l-ZI > Jxyj where J^y ~ 2/3J_l and for the XY case 
the sign of J is irrelevant the ground state was shown 
to be a vortex-liquid Mott insulator, a state with gapped 
magnetic excitation (excitations changing total ), how¬ 
ever having a gapless mode to non-magnetic excitation. 
For I J| < J!^y the ground state was shown to be a fully 
gapped non-degenerate state, adiabatically connected to 
the J = 0 case where the ground state is of simple prod¬ 
uct form and is composed of alternating = 0 com¬ 
ponents of triplets and singlets from rung to rung. The 
presence of multiple ground states when changing a single 
parameter (here ratio of exchanges along legs and rungs) 
is an indicator of frustration present in the system. 


STRONG RUNG-COUPLING LIMIT 


Let us start from the limit |J| <C J_l. For J = 0, 
even rungs in the ground state form rung-singlet states, 
whereas odd rungs form rung-triplet states. Hence for 
J = 0 the ground state manifold is 3^/^ times degenerate. 
Treating J perturbatively, integrating out singlets that 
occupy even rungs, we obtain an effective Hamiltonian 
describing a collection of odd rungs, 


3 L/2 

Tjl — 

- 1^ 


^T,T,+i+0(j4), 


( 2 ) 


WEAK RUNG-COUPLING LIMIT, 
BOSONIZATION 


In the other limit J ^ J± we use bosonization 
approach Q. We represent spin operators with the help 
of bosonic operators: 

S/j ^ a[Ji,L{x) + + (-l)^N,(x)]. (3) 


Decoupled chains have Neel-like quasi long-range or¬ 
der and the above representation captures the impor¬ 
tant low energy fluctuations by smooth bosonic fields, 
at wave-vector 0 and tt. Uniform spin magnetiza¬ 
tion is represented in terms of chiral currents 3i^l/r 
of a level-1 SU{2) Wess-Zumino-Witten model per¬ 
turbed by marginally irrelevant current-current interac¬ 
tions, describing an isolated antiferromagnetic Heisen¬ 
berg chain [l3|. 

We will need the following important operator product 
expansion (OPE) rules [l^. 


Ji,L/Ri^x)Ni{x' ,t') 


±i8abei{x', t') + i€abcNf{x', t') 
47r[u(T — t') ± i{x — x')] ’ 

(4) 


where v = 7rJ/2 is the spin-wave velocity of the 
Heisenberg spin-1/2 chain, known from the Bethe ansatz 
solution, and on the right hand side the dimerization 
operator ej, that is the continuum limit of (—l)'^SjjS;j_|_i 
has appeared. 

Treating the inter-chain coupling J±^ perturbatively for 
Jx! J ^ 1 in the continuum limit the staggered inter¬ 
chain coupling, H± = has the 

following form in terms of smooth bosonic fields. 


Hx 


J dx'Hx{x) 


Jx 


da:(Ji,L(a;) -f Ji,_R(a:))N2(a:). 


where Tj = Si_ 2 i-i-i + S 2 , 2 j-i-i are effective ^ = 1 spins 
formed along the odd rungs. Because the exchanges along 
both legs have equal strengths the lowest, second order 
in J contribution to the effective spin-1 chain formed 
on odd rungs vanishes. The ground state of our model 
Eq.JT]) for |J| <C Jx will be hence a direct product of 
singlets formed on even rungs and the ground state of a 
Heisenberg spin-1 chain formed on odd rungs. Depend¬ 
ing on the sign of J the ground state of this Heisenberg 
spin-1 chain is either ferromagnetic state (J < 0), or 
Haldane state [13 (J > 0). Eor the original Hamiltonian 
Eq.(IT]) ferromagnetic state of the effective S' = 1 chain 
Eq.(I2]) is the half-ferromagnetic state with the ground 
state total spin equal to half of the maximum possible 
value, S^ = L/2. This state will be called half-Eerro 
state. The effect of anisotropies on the Haldane state 
realized for 0 < J <C T_l has been studied recently [111- 
Eor the rung-alternated ladder with J > 0 it was shown 
that the application of an external magnetic field induces 
the half-magnetization plateau state 121. 


+Jx J d:E(J 2 .L(a:)+J 2 .fl(a^))Ni(x). (5) 

The Hx perturbation has non-zero conformal spin and 
does not open a gap in the first order of J_l . The relevant 
scalar operator from inter-chain exchange comes in the 
second order of Jj_ coupling. 

Using OPE for the same-leg operators at short- 
distances Eq.Q and integrating with the relative coor¬ 
dinates we obtain in the order the following relevant 
contributions that should be added to effective Hamilto¬ 
nian describing the long wavelength properties of decou¬ 
pled bosonic chains, 

~ - Ji y dx[3ei(a;)e2(a;) - 2Ni(a;)N2(x)]. (6) 

At this point it is convenient to introduce 4 Majorana 
fermions [l^ . The perturbation that we identified in 
Eq.® is translated as mass term of triplet and singlet 
Majorana fermions and the complete ladder Hamiltonian 
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Half-Ferro | Hi ? Rung-Singlet 
J <0 0 J«Ji J»J^ J 

FIG. 2. Analytically conjectured ground-state phase di¬ 
agram of rung-alternated Heisenberg ladder as function of 
J. Hi denotes a Haldane phase, however only half of the 
rungs (odd numbered rungs) provide the effective S = 1 spins, 
whereas the other half (even numbered rungs) are in the ap¬ 
proximate rung-singlet states and mediate antiferromagnetic 
exchange among effective S' = 1 spins. 
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FIG. 3. Geometry of honeycomb-ladder model. For a = 0 
rung-alternated ladder presenting in Fig. 1 is recovered. 


HONEYCOMB-LADDER MODEL 


in Majorana basis looks, 


Hm= 


/' 




+ I (7) 


where the masses of Majoranas are: mt ^ 5J'j_ and 
nis ^ 3J]_. Following the analyses of the ground states of 
the original ladder system from Majoranas [IQ we can 
conclude that the ladder system for J^/J —)► 0 is in the 
phase that is adiabatically connected to the rung-singlet 
phase of the uniform ladder with antiferromagnetic ex¬ 
changes. This is surprising, because the spin exchanges 
are ferromagnetic along every other rung. 

Comparing the effective Hamiltonian Eq. (ED and the 
original one Eq. (HD reveals a very important compe¬ 
tition forJ > 0. Namely, for small 0 < J <C J± the 
next nearest neighbour spins on the same chain belong¬ 
ing to the odd rungs show antiferomagnatictendencies, 
{^i, 2 j+i^i, 2 j+ 3 ) < 0, whereas for the single-chain domi¬ 
nated regime J » J_l, it is clear that {SijSij+ 2 ) > 0 for 
any I and j. There is no such competition on J < 0 side 
and hence the case of ferromagnetic legs is simpler and 
for the entire region of J < 0 the ferromagnetic phase of 
the effective spin-1 chain is expected to be realized (for 
the original ladder it is the half-Ferro state with ground 
state total spin = L/2). 


In this section we will study a slightly generalized case 
of the rung-alternated ladder, introducing the additional 
parameter a that relaxes the condition of equal absolute 
values of exchanges along even and odd rungs. This way 
we obtain a SU{2) symmetric version of the so called 
honeycomb-ladder model [l6| . 

L L 

H = j^ -( 8 ) 

i=j j=l,l 

To simplify notations in the following we will set J± = 1. 

For a G (—1,1) the spin system is frustrated, because 
the number of antiferromagnetic bonds per ladder pla- 
quette is odd. We will study the ground-state phases for 
arbitary values of a and J (both negative and positive). 
The simpler case to start from is J <0. 


Honeycomb-ladder with ferromagnetic legs 


For J < 0 we can estimate the boundary of the fer¬ 
romagnetic phase. We derive the single-particle disper¬ 
sions, composed of 4 branches, since the unit cell con¬ 
tains 4 spins. 


ei(k) 

£3(k) 


S4(k) 


— -F J cos fc, £ 2 (fc) = — — J cos k 
a \J2 + J'^ + cos2fc 

2 V 2 

a ^ y/2 + -\- J'^ cos 2k 

"2 + 7! ■ 


(9) 


So far, using analytical approaches, we have estab¬ 
lished ground-state phases of rung-alternated Heisenberg 
spin-1/2 ladder in the limiting cases of: 0 < J J_l and 
J » Jj_, where a Haldane phase of effective S' = 1 spins 
formed along the odd rungs and a rung-singlet phases 
are stabilized respectively. The conjectured sequence of 
phases with changing J is depicted in Fig. [2j 

To check whether there are additional phases for in¬ 
termediate values of J ^ Jj_ or there is a direct phase 
transition between Hi and RS states, we will use numer¬ 
ical approaches. 


the chemical potential is attached to the minimum of 
£i(fc) band, realized at fc = 0, /r = a/2 + J■ By looking 
at the single magnon instability we estimate the bound¬ 
ary of the fully polarized state, equating the chemical 
potential to the minimum of e^^k) band, which is given 
for any system size by 


OtpM = — 




'2\J\ 


for I J| » 1. 


( 10 ) 


Next, for | J| 1 the single-particle dispersion bands 

become flat and we can estimate the transition from the 
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ferrimagnetic state into the half-Ferro state with total 
S' = L/2 by equating the chemical potential with the 
maximum of £ 3 (fc) band realized at fc = ±7r/2, 

ac = - J -1 + 0{J‘^). (11) 

Transitions induced by changing a < 0 can be under¬ 
stood in a simple way for |J| <C 1. In the half-Ferro 
state spins on even sites form rung-singlets, whereas on 
odd sites spins form effective spins S = 1. The effective 
spin-1 chain, formed by spins on odd sites, is in the fully 
polarized phase as follows for J < 0 from Eq. [2j The in¬ 
teraction between the spins belonging to even sites along 
the leg direction is mediated by intermediate S' = 1 spins 
(all of which are pouting in the same, spontaneously cho¬ 
sen, direction due to the fact that they are in fully po¬ 
larized ground state) and this interaction is ferro since 
J < 0. Hence, the system is equivalent to the direct prod¬ 
uct of the ferromagnetic state of the spin-1 chain (formed 
by spins belonging to odd sites) and a 2-leg spin-ladder 
with antiferromagnetic exhcange along rung and ferro¬ 
magnetic exchange along legs[17[ (the ladder is formed 
by spins belonging to even sites). Decreasing a weak¬ 
ens the antiferromagnetic coupling along the rung and 
the effective magnetic field (produced by the fully polar¬ 
ized neighbouring spin-1 chain) induces two consequtive 
second order commensurate-incommensurate transitions, 
first one from rung-singlet state of the two-leg ladder 
(with ferromagnetic legs and antiferromagnetic rungs) to 
an intermediate Luttinger liquid state with finite polar¬ 
ization and then to fully polarized state [l8j|. 

Single-magnon (equivalently hard-core boson) disper¬ 
sions as presented in Fig. |4]shed light on the evolution 
of the low-energy excitation spectrum of model Eq.® 
as function of a for a < 0 and J < 0 and fully con¬ 
firm the above picture anticipated from the interpene¬ 
trating spin-1 chain and the two-leg ladder with ferro¬ 
magnetic legs and antiferromagnetic rung exchanges. For 
a < — vT+T7^— J system is in the ferromagnetic state 
and the low energy excitation is a conventional fer¬ 
romagnetic magnon, with ^ cos/c-like gapless quadratic 
dispersion shown in Fig. |4] (a). At the boundary be¬ 
tween ferromagnetic and ferrimagnetic phases there are 
two quadratic dispersions, two kinds of magnons with 
the mass ratio of y/2 presented in Fig. |4](b). Inside the 
ferrimagnetic state the ground state total spin changes 
continuously L/2 < S’^ < L and there are two differ¬ 
ent kinds of dispersions: one quadratic at low momenta 
and a second one that is linear shown at zLkp it Fig. 0] 
(c), giving Luttinger liquid like properties. At the bound¬ 
ary between ferrimagnetic and half-Ferro states again two 
quadratic in momenta low energy dispersions are present 
shown in Fig. 0] (d) and inside the half-Ferro state with 
S'^ = L/2 only one magnon branch with ^ cos/c-like 
quadratic low energy dispersion remains as shown in Fig. 
0](e). Both phase transitions, from ferromagnetic to ferri¬ 
magnetic and from ferrimagnetic to half-Ferro states are 


e(k) _e(k) 





(b) 








-71/2 0 k 71/2 -71/2 0 k 7 c/2 
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FIG. 4. (Color online) Single particle dispersion relations 
(a) deep in ferro state, (b) at the boundary of ferro and ferri 
states, (c) inside ferri state, (d) at the boundary of ferri and 
half-ferro states and (e) inside the half-Ferro state. the 
dashed line indicates the chemical potential. 

second order commensurate-incommensurate phase tran¬ 
sitions, where the linear mode disappears in favor of a 
quadratic dispersion. The overall gapless quadratic mode 
remains in the background, due to the spontaneously 
broken SU(2) symmetry in all phases where the ground 
state is not a global spin singlet. This makes SU{2) sym¬ 
metric honeycomb-ladder model Eq.® for J < 0 a very 
attractive and simple case to study the behavior known 
as unsaturated ferromagnetism, an effect that has been 
noticed to occur in other frustrated systems [l9l - E^ . With 
further increase of a to positive values there is a phase 
transition from half-Eerro to rung-singlet state, not cap¬ 
tured by the single particle picture. 


Effective model for a ~ 1 and | J| -C 1 

Eor a ~ I and |J| <C 1 we can derive an effective 
model for spins on odd rungs by integrating out spins 
belonging to even rungs that are in rung-singlet states. 
To capture the difference between J < 0 and J > 0 cases 
we have to go beyond the lowest (second) order in inter¬ 
rung coupling J that we treat in perturbation theory. To 
the third order in J the effective ladder model formed by 
spins belonging to odd rungs is given by the the following 
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FIG. 5. Effective model for spins belonging to odd rungs 
valid for a ~ 1 and J 1, where spins belonging to even 
rungs form approximate rung-singlet states. Phase transitions 
are expected for Jr < 0. 


hamiltonian (see Fig. [5]), 


L/2 

Heff = J\\ ^[Sl,2i-lSl,2i-|-l + S2.2j-lS2,2j-|-l] 

J = 1 
LI2 

+ ^^[Sl,2j-lS2,2j-|-l + Si,2j-|-lS2,2i-l] 

i=i 

L/2 

+ Sl,2j-t-lS2,2/-|-l (12) 

i=i 


where 


2(a-bl) 4(Q;-bl)2 

J2 3J3 

2(Q;-b 1) 4(Q;-b 1)2 

a2 J2 - 1 3 J3 

a + 1 2{a -I- 1)2 


(13) 


Interchanging legs with diagonals of the effective 
model presented in Fig. [S]by interchanging spins on ev¬ 
ery other rung, the model Eg. (1121) for parameters given 
in Eq. (1131) is equivalent to a two-leg spin ladder with 
antiferromagnetic legs, ferromagnetic diagonals and rung 
exchange that changes from antiferromagnetic to ferro¬ 
magnetic with decreasing a. For the two-leg antiferro¬ 
magnetic ladder weakly coupled by competing diagonal 
and rung exchanges (where bosonization is applicable) 
three phases are expected to be stabilized with decreasing 
a : rung-singlet, dimer and Haldane phase [l^ . Parame¬ 
ters of our effective model are outside the weak-coupling 
limit, but later with the help of numerical simulations 
we will show that the same sequence of phases are also 
realized in our effective model Eq.([T2|). In particular, the 
dimerization pattern of the original ladder model will be 
dimers formed along next nearest neighbor diagonals, in¬ 
volving spins belonging to odd rungs. 


Vicinity ofa~—1, 0<J‘Cl 


For J > 0 and 0 < J <C 1 we can as well estimate pos¬ 
sible ground states around a ~ — 1, where we can borrow 
the results from the mixed diamond chain [illiil. For 
J = 0 and a = — 1 the spins belonging to even rungs 
are disconnected and spins belonging to odd rungs in the 
ground state are in the spin-triplet configuration, form¬ 
ing S' = 1 spins. For 0 < J -C 1 there is a competition in 
the nature of the exchange between the spins belonging 
to the even rungs: for a > —1 the direct exchange S = 1 
is antiferromagnetic, whereas the exchange mediated by 
nearby S = 1 spins is ferromagnetic. For a < —1 there is 
no such competition and effective spins S = 1 are formed 
on each rung. 

One possibility that the above mentioned competition 
for a > — 1 gets resolved is that some of the even rungs 
choose to be in triplet state and others in singlet state pe¬ 
riodically alternating as happens in the mixed-diamond 
chain |24| where consecutive odd number of M rung- 
triplets (coupled antiferromagnetically with each other 
by J > 0) will be sandwiched between the rung-singlets. 
Coupling an odd number M oi S = \ spins by antiferro¬ 
magnetic exchange and assuming open boundary condi¬ 
tions, the M-rung segment will be in the triplet state in 
the ground state, forming an effective S = 1 spin. The 
approximate rung singlets in the case of rung-alternated 
ladder (as opposed to the exact rung singlets realized in 
mixed-diamond chain 2^ that cut the chain) will medi¬ 
ate an effective antiferromagnetic exchange among the 
above mentioned effective S' = 1 spins formed by M-rung 
segments, giving rise to generalized Haldane states with 
an enlarged unit cell composed of M -|- 1 ladder plaque- 
ttes. Such a state for M = 3 is depicted in Fig. Inland 
called Haldane-dimer. 

There are in total 5 different ground states in the 
mixed-diamond chain when changing the equaivalent of 
a from a > —1 to a < —1 with M = 1, 3, 5, 7 and M = 
oo. The M = 1 case is equivalent to the Hi state, 
a Haldane state of S = 1 spins formed on odd rungs, 
and the M = oo case is equivalent to the Haldane state 
of the effective S = 1 spins formed on every rung. For 
the spin-1/2 ladder with alternated rungs M = 1,3,5 
and M = oo are suggested to have hnite extent in the 
presence of exchange anisotropy |25l| . 


NUMERICAL RESULTS 


Rung-alternated ladder 


Initially we will present our numerical data for a = 0 
case corresponding to the rung-alternated ladder model 
Eq. (1). We use both Lanczos simulation and the density 
matrix renormalization group (DMRG) approach 261.1271 
in order to access large system sizes. For the case of 
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FIG. 6. Cartoon of one of the possible ground state config¬ 
urations of the Haldane-Dimer phase that can be realized 
for a~—1,0<J<C1. Spins encircled by open rectangles 
form approximate rung-singlet state, whereas those encircled 
by shaded rectangles form rung-triplet states. 6 spins encir¬ 
cled by the dotted rectangle form effective S = 1 spins which 
are connected via intermediate singlets to produce an effec¬ 
tive Haldane chain. In the mixed-diamond chain the singlets 
depicted above are exact eigenstates and they do not mediate 
any exchange among the effective S = 1 spins. 


ferromagnetic legs we systematically obtain (both using 
large scale DMRG as well as Lanczos algorithm for both 
periodic and open boundary conditions) that the ground 
state belongs to the multiplet with total spin half of the 
maximal possible value, S’^ = L/2 for any J < 0. 

In the following we will discuss J > 0, where we pre¬ 
dicted at least two different phases in the limiting cases 
J <C 1 and J ;g> 1 respectively. In Fig. [7] we plot the fi¬ 
delity susceptibility 2^ ^ | with changing control param¬ 
eter J for different system sizes. 


2 \n\{MJ)\MJ + ^J))\ 

L 5 ™0 (5J)2 


(14) 


where |'0o(<7)) is the (non-degenerate) ground state wave- 
function for the corresponding parameter J. We see that 
there is a well pronounced peak in the fidelity susceptibil¬ 
ity and the height of the peak increases with system size, 
whereas the width decreases. We extrapolate the loca¬ 
tion ot the peak to J = Jd ~ 0.45 in the thermodynamic 
limit. Thus, we can estimate the extent of the Hi phase 
for a = 0as0<J< Jd. Note, the rather similar esti¬ 
mate of Jd follows from the position of the level-crossing 
of the lowest excited states, which are triplet with mo¬ 
mentum fc = TT in Hi phase (J < Jd) and triplet with 
/c = 0 in its neighboring phase (J > Jd). 


Honeycomb-ladder with ferromagnetic legs: J < 0 

We start by presenting our numerical data with the 
case of ferromagnetic legs. In Fig. |8] we plot for the 
honeycomb-ladder with ferromagnetic legs the ground 
state multiplicity as function of a for J = — 1, which is a 
typical behavior in the whole J < 0 region. The bound¬ 
ary of the ferromagnetic phase is captured exactly from 
the spin-wave instability (indicated by dashed vertical 



Q nfiwnrannnl a l l CiOOPm*'' ^ i _|_|_ i n P P ■ 

0.3 0.35 0.4 0.45 ^ 0.5 0.55 0.6 


FIG. 7. (Golor online) Fidelity susceptibility of the rung- 
alternated ladder with antiferromagnetic legs for periodic 
boundary conditions and 3 different system sizes. In DMRG 
calculations periodic boundary conditions are assumed (re¬ 
stricting considerably the available ladder lengths) to avoid 
degeneracies of the Haldane-like states due to edge spins for 
the open boundaries. 



FIG. 8. (Golor online) DMRG results for the lowest energy 
levels in different total subspaces relative to the lowest 
energy in the — 0 subspace as function of a for J — —1 
and L = 96 rungs for open boundary conditions. The dashed 
line indicates the boundary of the ferromagnetic state. 


line in Fig. |8]). Here we only present the data from which 
we determine the boundaries of the half-Ferro state. 

We see from this plot that for J = — 1 the half-Ferro 
state S'^ = L/2 is sandwiched between a ~ —0.2 and 
a ~ 0.36. Note, for 0.36 < a < 0.4 the excitations from 
the singlet ground state to low total spin states e.g. to 
states with S'^ = 1 and 2 become practically gapless, thus 
we cannot rule out the existence of an intermediate thin 
phase between half-Ferro and rung-singlet states based 
on our numerical data. 

For values of a < —0.2 energies of the lowest states 
with > L/2 merge gradually with the ground state 
(only one state = L/2 + 1 is indicated in Fig. [5]) 
until the energy of the fully polarized state with = 
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0.9 0.92 0.94 0.96 0.98 a 

FIG. 9. (Color online) Lanczos results of the lowest excited 
states relative to the ground state obtained for L = 12 rungs 
of the effective model Eg. (1121) for (a) J = 0.1 and (b) J=0.2. 
A similar picture is expected to hold for the Honeycomb- 
ladder model with L = 24 rungs for the same extent of a 
and J. Periodic boundary conditions are used that allows to 
assign a definite lattice momentum to each level. 


L becomes degenerate with the ground state energy for 
a < apM, ocfm{J = ~1) = 1 ~ — —0.414. 


Honeycomb-ladder with antiferromagnetic legs: 

J > 0 

For J > 0 case we start presenting numerical data near 
the point a ~ 1 for small J with changing a. To distin¬ 
guish different phases it is usefull to start from looking 
at the gap between the ground state and lowest excited 
states. In Fig. Owe depict the lowest excited states as 
function of a. We use the effective model Eq. (HU to 
reach system sizes of L = 12 rungs, that is equivalent 
to L = 24 rungs for the Honeycomb-ladder model Eq. 
(|8|). We have checked that for available system size (up 
to L = 12 for Honeycomb-ladder model) agreement be¬ 
tween the low energy levels of effective and full models is 
perfect for small J values. In FigOwe present the level 
spectroscopy results for J = 0.1 (a) and J = 0.2 (b) for 
the effective model Eq. ([8]). In an antiferromagnetic lad¬ 
der with a uniform antiferromagnetic exchange the lowest 
excited state in the rung-singlet phase is a triplet state 
with wavevector /c = tt in units of the ladder lattice con¬ 


stant. Since in our model the unit cell is made of two pla- 
quettes, in the rung-singlet phase of Honeycomb-ladder 
(in the phase that is adiabatically connected with rung- 
singlet phase of the uniform ladder, but with a unit cell 
half of the Honeycomb-ladder model) the lowest excited 
triplet should have momentum /c = 27r in the units of the 
Honeycomb-ladder unit cell that is equivalent to fc = 0 
momentum. 

We see that with decreasing a below a = I the gap to 
the lowest excitation (triplet state with fc = 0 momen¬ 
tum) shows a minimum and then with reducing a this 
lowest triplet excitation level crosses with the lowest ex¬ 
cited singlet state that has momentum fc = tt. Note, for 
any J > 0 and any a the ground state is a spin singlet 
state with fc = 0 momentum. There is a finite extent in 
a where the lowest excited state is a singlet state with 
fc = TT. With further reducing a there is a level-crossing 
between the lowest spin singlet excitation with fc = tt 
and spin triplet excitation with fc = tt. The spin triplet 
excitation with fc = tt in units of Honeycomb-ladder unit 
cell is the lowest excitation on top of the Haldane state 
that is defined on the effective spin-1 chain with the same 
unit cell as the original microscopic model. One can use 
the abovementioned two level crossings in excited states 
to estimate the stability region of the intermediate dimer 
phase. In fact with increasing system size the singlet ex¬ 
citation at fc = TT should get degenerate with the ground 
state singlet in the dimerized phase. The boundary be¬ 
tween the rung-singlet and dimer states can be estimated 
from the position of the minimum of the gap of the fc = 0 
triplet state. 

Comparing the energy levels of the effective ladder 
model Eq. (113 for different system sizes with L < 12 
rungs we see that the energy of singlet state with fc = tt 
momentum decreases faster, with increasing the system 
size, than energies of the triplet states in the parameter 
region where the dimerized phase is expected. 

The numerical ground-state phase diagram of the 
honeycomb-ladder model obtained with the help of 
DMRG simulations is presented in Fig. [TOl For J < 0 
there are four different phases realized with decreasing a: 
rung-singlet, half-ferro, ferrimagnetic and ferromagnetic. 
For the case J > 0 the rung-singlet state, the Hi and 
the conventional Haldane state and different dimerized 
states: NNND and Haldane-Dimer are realized. Between 
the NNND-Dimer and Haldane-Dimer states we can not 
locate numerically the phase transition line, neither can 
we exclude emergence of an intermediate (gapless) state 
located around a = 0. 

It is worth noting that the topology of the Hi phase 
realized for J > 0 can be captured by studying one pla- 
quette of the ladder, L = 2. Consider e.g. the a = 0 
case. For this case for J < 1/^/2 a triplet state is realized 
as ground state, whereas for J > l/\/2 the ground state 
becomes a singlet, a direct product of the singlet states 
on the first chain (2-site chain) and on the second chain 




















FIG. 10. Numerical ground-state phase diagram of 

the honeycomb-ladder model in the parameter plane 
NNND-Dimer stands for next nearest neighbor diagonal dimer 
phase where dimers are formed along next nearest neighbor 
diagonals involving spins of odd rungs. In the vicinity of 
a = — 1 and for 0 < J <C 1 before the transition from 
Haldane-Dimer to Haldane phase additional phases may occur 
(e.g. with M = 5 and M = 7 as discussed in previous sec¬ 
tions). In dimer phases ground states are doubly degenerate 
in the thermodynamic 

limit. 

(that is an exact eigenstate for any J in the case of a 
single plaquette). Hence at J = l/-\/2 there is a triplet- 
singlet level crossing in the ground state of one plaquette. 
The (threefold) degeneracy of the ground state for small 
values of J is a particular case and omit:it stems from 
the fact that there is only one effective spin 1 (formed 
on one of the two rungs). As soon as the number of 
ladder plaquettes is increased and more than one effec¬ 
tive spin 1 is formed on odd rungs the ground state 
becomes a singlet (for periodic boundary conditions) for 
the whole range of J > 0 and there is no level crossing in 
the ground state any more. However, when one assem¬ 
bles many plaquettes into the ladder geometry, instead 
of the level-crossing, one can identify the avoided level 
crossing in the lowest energy singlet states of the finite 
ladder (for system sizes L < 12 rungs), that is located at 
J ^ 0.5 (data not shown). 

The phase transition points indicated in Fig. [in] for 
J > 0 were obtained by studying the behavior of the fi¬ 
delity susceptibility as function of a for different values of 
J as presented in Fig. [TT] For small values of J (roughly 
J < 0.5) the fidelity susceptibility shows typically well 
pronounced four peaks, whereas for J > 0.5 only two 
peaks are visible, one for positive and a second one for 
negative a. The peak for the a > 0 side becomes less 
and less pronounced with increasing J > 1. 

To describe the regime corresponding to J 1 in 
Fig HD we present the behavior of the lowest excitation 
gap as function of a using DMRG for large value of J = 5. 
In order to access large system sizes we use open bound¬ 
ary conditions. One can see that with decreasing a first 



-0.5 0 CK 0.5 1 



FIG. 11. (Color online) Ground state fidelity susceptibility 
per site as function of a for (a) J = 0.4 and for J = 1 showing 
two peaks. In DMRG calculations periodic boundary condi¬ 
tions are assumed to avoid degeneracies of the Haldane-like 
states due to edge spins for open boundaries. 

there is a minimum in gap and then there is a cusp¬ 
like behavior. Using the finite system size data for sys¬ 
tems with L = 48,96 and 144 rungs the position of the 
gap minimum in the thermodynamic limit extrapolates 
clearly to negative values of a. Starting from the rung- 
singlet phase, the gap decreases linearly with decreasing 
a and the position of the minimum of the gap we inter¬ 
pret as a boundary of the rung-singlet phase. 

On the other hand, for J ^ 1, extending the bosoniza- 
tion analyses to a ^ 0 gives that for a > 0 the 
rung-singlet phase smoothly evolves into the rung-singlet 
phase of the uniform antiferromagnetic ladder realized for 
a ^ 1. For Of < 0 interestingly bosonization suggests the 
sequence of two consecutive second order phase transi¬ 
tions, first from rung-singlet to an intermediate dimer 
phase and then from dimer to Haldane phase with de¬ 
creasing a. Hence we expect to see two values of a < 0 
where gap should close in the thermodynamic limit. In¬ 
stead we see only one minimum in the finite-size gap data 
presented in Fig. [121 The reason why we do not see the 
second minimum may be the fact that finite-size effects 
are still large (even for L = 144 rungs). In addition, 
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FIG. 12. Gap between the ground state and the first excited 
(triplet) state as function of a for J = 5 obtained by DMRG 
using open boundary conditions. 

since we use open boundary conditions, it is difficult to 
separate true bulk gap from the boundary gap of the 
Haldane phase (realized on the left side from the kink 
in Fig. [T^ . It is desirable to study the gap for peri¬ 
odic boundary conditions, however DMRG calculations 
become less accurate and only much smaller system sizes 
can be addressed. 


CONCLUSIONS 

We have studied the ground-state phase diagram of 
the rung-alternated SU{2) symmetric spin—1/2 ladder. 
Both cases with ferromagnetic as well as antiferromag¬ 
netic leg exchanges have been considered. For the case of 
ferromagnetic legs we showed that a unique ferrimagnetic 
ground state emerges, with ground state magnetization 
equal to half of the maximum possible value, for arbitrary 
strength of the leg exchanges. The case of antiferromag¬ 
netic leg exchange is much richer and depending on the 
ratio of leg to rung couplings several different ground 
states can emerge starting from the Haldane phase Hi 
for small leg couplings and ending with the rung-singlet 
phase for strong leg couplings. Based on Fig. (TU] it is 
tempting to speculate that dimer order extends to a = 0 
and hence the intermediate phase of rung-alternated lad¬ 
der can be dimerized, even though we have not succeeded 
in either directly measuring dimerization order, or find¬ 
ing a second singlet state as the lowest excited state of 
the finite chain or even resolving a finite excitation gap 
numerically. 

We have as well studied a generalization of rung- 
alternated ladder: the spin—1/2 Heisenberg system on 
honeycomb-ladder lattice. For the case of ferromag¬ 
netic legs we have identified a peculiar Luttinger liquid 


ferrimagnetic state, where the ground state magnetiza¬ 
tion changes continuously as function of system param¬ 
eters and low energy gapless excitations consist of two 
branches one of which is linear and another quadratic 
in momentum. For the case of antiferromagnetic leg 
couplings different short-range ground states, including 
those with possible Haldane-like topological order have 
been suggested to occur. 
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